Magnetism, structure, and charge correlation at a pressure-induced Mott-Hubbard 

insulator-metal transition 



Yejun Feng, 1,2 R. Jaramillo, 3 A. Banerjee, 2 J. M. Honig, 4 and T. F. Rosenbaum 2 

The Advanced Photon Source, Argonne National Laboratory, Argonne, IL 60439, USA 
2 The James Franck Institute and Department of Physics, 
The University of Chicago, Chicago, IL 60637, USA 
'School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 01238, USA 
* Department of Chemistry, Purdue University, West Lafayette, IN 47907, USA 

We use synchrotron x-ray diffraction and electrical transport under pressure to probe both the 
magnetism and the structure of single crystal N1S2 across its Mott-Hubbard transition. In the 
insulator, the low-temperature antiferromagnetic order results from superexchange among correlated 
electrons and couples to a (1/2, 1/2, 1/2) superlattice distortion. Applying pressure suppresses the 
insulating state, but enhances the magnetism as the superexchange increases with decreasing lattice 
constant. By comparing our results under pressure to previous studies of doped crystals we show 
that this dependence of the magnetism on the lattice constant is consistent for both band broadening 
and band filling. In the high pressure metallic phase the lattice symmetry is reduced from cubic to 
monoclinic, pointing to the primary influence of charge correlations at the transition. There exists a 
wide regime of phase separation that may be a general characteristic of correlated quantum matter. 

PACS numbers: 71.30.+h, 75.30.Kz, 61.05.cp, 64.75. Qr 



I. INTRODUCTION 

One of the great challenges of understanding corre- 
lated materials is teasing apart the relative influences 
of the spin, charge, orbital, and lattice degrees of free- 
dom. The point where an insulator becomes a metal 
highlights acutely the competition between mechanisms, 
but it also affords a special opportunity to limn the perti- 
nent physics when different routes across the phase tran- 
sition are available. Going back to the original ideas 
of Mott and Hubbard, we know that strong Coulomb 
repulsion between electrons on a single lattice site can 
localize ch arge even when band theory predicts metal- 
lic behavioiQEI, anc } antiferromagnetism was attributed 
to a consequence of superexchange between localized 
electrons^. At the same time, Slater claimed that an- 
tiferromagnetism alone could account for the formation 
of the insulating gapR All can be subsumed by symmetry 
changes wrought by a structural phase transition^. 

The cubic pyrite crystal MS2 has long been reco gniz ed 
as a canonical Mott-Hubbard correlated insulatop^EHISl 
Band structure calculations^ put the sulfur 3p band and 
Ni tig band well below the half filled Ni e g band, pointing 
to on-site Coulomb repulsion as the source of the insu- 
lating energy gap E g , which lies in the range 1-10 meV 
(Ref. |12|). The small size of this gap demonstrates that 
N1S2 is an incipient Mott insulator^ with the Coulomb 
repulsion U comparable to the e g bandwidth, W = 2.1 
eV (R ef. r[T] b The gap can be suppressed either by Se 
^ p m ^2i8ii4j or a ppii ec i pressur^UHlB k u t the scale of the 

pressure required to drive the gap to zero in the pure limit 
has introduced technical obstacles to systematic stud- 
ies of the competition between electronic, magnetic, and 
structural correlations at the quantum phase transition. 
Doping with Se expands the lattice and reduces the Ni 3c? 
bandwidth^!, and the transition in this case is thought 



to be driven by increasing charge transfer between the 
Ni 3g! and Se 4p bands^. Applying pressure tunes the 
ratio U/W and provides a more direct approach to the 
Mott-Hubbard model. 

We use synchrotron x-ray diffraction and electrical 
transport in a diamond anvil cell to parse the roles of 
the low temperature antiferromagnetism and the lattice 
structure, both through the insulating state and at the 
transition to the metal in the Mott-Hubbard system, 
N1S2. The reduced symmetry in the metal to monoclinic 
- a highly unusual occurrence for correlated materials de- 
scribed below - eliminates the change in the structure as 
a likely origin of the derealization of charge. By compar- 
ison of the pressure-induced transition in the pure com- 
pound to the insulator-metal transition driven by chem- 
ical substitution of Se for S, we identify the charge de- 
grees of freedom as the predominant driving mechanism. 
Realizing the insulator-metal transition in high quality 
single crystals of a stoichiometric material using applied 
pressure further clarifies the physics by avoiding compli- 
cations that arise from chemical disorder, most notably 
the competition between Anderson localization and the 
Mott transition. 



II. EXPERIMENTAL METHODS 

High-pressure single-crystal x-ray diffraction measure- 
ments were carried out at bcamlincs 4-ID-D and 6-ID-B 
of the Advanced Photon Source. In a vertical scattering 
geometry with a psi-diffractometer, a high g-resolution 
(FWHM - 1 x 1CT 3 A~ ) is achieved using a 50 [im 
size detector slits positioned 1.3 m away from the sample 
along the 28 arrrP. The use of double-bounce Pd mirrors 
for 20 keV x-rays and an energy discriminating Nal scin- 
tillation detector eliminated higher-harmonic contamina- 
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tion of the diffraction signal. Our single crystals were 
grown by the Te flux method to remove potential com- 
plications from excess impurity concentrations^. Crystals 
were 25 to 50 /im in diameter and fit well within the dia- 
mond anvil cell pressure chamber. Five different crystals 
were studied under pressure using a methanohethanol 4:1 
mixture for the pressure medium. Base temperature var- 
ied between 3.5 and 5.8 K, and the pressure was cali- 
brated in situ using silver diffractiorP. 

The problem of vacancies, common to sulfides, is well 
characterized^ in N1S2. Both S (about 4%) and Ni (vary- 
ing) vacancies can be determined from the measured lat- 
tice constants and electrical resistivity; in this way we 
estimate our sample stoichiometry to be NiSi.96- There 
are two coexisting antiferromagnetic structures in NiS2- 
The Ml antiferromagnetic order with a wavevector (1, 0, 
0) emerges from a second order phase transition at Tjvi = 
37 to 54 K, where T/vi strongly depends on vacancy con- 
centration and varies from sample to sample (i.e. Ref. 
[PJ). The M2 antiferromagnetic order has a (1/2, 1/2, 
1/2) wavevector and emerges at a first order transition 
at Tjv2 = 30 ± 1 K, where the transition temperature 
is consistent across all published reports including those 
for different vacancy concentrations. However, vacan- 
cies are responsible for the variable canting angle in an 
antiferromagnelP^, as observed in the M2 phase of N1S2 
(Ref. [SHU]). 



III. ANTIFERROMAGNETISM AT AMBIENT 
AND HIGH PRESSURES 



One outstanding question in the field of Mott- Hubbard 
systems concerns the role of period doubling antifer- 
romagnetism at the insulator-metal transition. The 
magnetostrictive response of the lattice to M2 is likely 
rhombohedraP^ and the distortion from cubic symme- 
try is extremely small, with a relative lattice constant 
change Aa/a - 2 x 10~ 4 (Ref. [2D]). Therefore, for 
our diffraction measurements, we were able to model 
the insulating ground state using a cubic matrix. At 
ambient pressure and T < Tjv2 we report the discov- 
ery of charge-originated superlattice diffraction peaks at 
(1/2, 1/2, l/2)-type positions in reciprocal space (Fig. 
1), corresponding to the M2 magnetic structure. Given 
the compatibility of diamond anvil cell technology with 
high energy x-ray diffraction, this opens up the possibil- 
ity of combined magnetic and structural studies of the 
pressure-driven Mott-Hubbard transition in pure NiS2- 
The temperature dependence of the superlattice intensity 
(Fig. la) scales linearly with the magnetic M2 neutron 
diffraction intensity^. The x-ray supperlattice and neu- 
tron M2 diffraction intensities scale quadratically with 
the superlattice displacement and also the magnetic mo- 
ment, respectively. The scaling in Fig. la thus points to 
a linear coupling between the M2 magnetic moment and 
the superlattice displacement. A linear coupling between 
magnetism and the lattice is rather common in solids, 
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FIG. 1: (color online), (a) Superlattice intensity Isl at am- 
bient P plotted as a function of temperature. Isl is averaged 
over (3.5, 1.5, 1.5), (3.5, 0.5, 1.5), (2.5, 1.5, 1.5), and (2.5, 
0.5, 1.5), and normalized to (3, 1, 1). Data were taken while 
warming (pink) and cooling (blue) through Tjv2 = 29.2 K. 
The x-ray superlattice intensity scales linearly with the in- 
tensity of magnetic neutron diffraction from the (1/2, 1/2, 
1/2) M2 magnetic order (open circles, data from Ref. [9]), as 
well as the linear thermal expansion \AL/L\ along (1, 1, 1) 
(black solid squares, data from Ref. [20 ). Solid line is a guide 
to the eye. (inset) H, K, L scans of (1/2, 1/2, 1/2) superlat- 
tice diffraction at ambient P and T = 5.8 K. (b) Resistivity 
of NiS2 at several pressures. All data follow an Arrhenius 
form through the first order phase transition at Tjv2, which 
is marked by a kink, (c) P-x-T phase diagram of the M2 
magnetic phase. Tn2 is plotted as a function of low temper- 
ature cubic lattice constant ao for both MS2 under pressure 
and NiS2- ;E Se ;E (blue open symbols: [14.16 17]; blue solid cir- 
cle: our transport measurements). The two dot-dash lines 
mark the first-order boundaries of the insulator-metal transi- 
tion driven by pressure and chemical substitution, and bound 
the M2 order. 



and most often observed in low-dimensional syste ms^. 
By comparison, the external magnetostrictiorpSlSSJ, mea- 
sured by the linear thermal expansion, scales quadrati- 
cally with the M2 magnetic moment (Fig. la). 

The intensity of the superlattice peaks provides a quan- 
titative estimate of lattice distortion d through the rela- 
tion I SL / I(3ii) = \q-S/2\ 2 S(SL)/S(31l), where q is the 
superlattice ordering wavevector and S(SL) and 5(311) 
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FIG. 2: (color online), (a) Comparison of the high-resolution 
longitudinal linescans of the superlattice at four pressures be- 
tween and 2.84 GPa. Data are plotted from the center qo of 
each individual peak and displaced vertically with scale bars 
representing individual intensity relative to that of the (3, 
1, 1) order. All lineshapes are nearly resolution limited ap- 
proaching the phase boundary. This indicates that the crystal 
remains cubic with long M2 correlation lengths throughout 
the insulating phase. The antiferromagnetic M2 order dis- 
appears at higher pressures as exemplified by two null scans 
through the superlattice positions at 2.90±0.05 and 3.86±0.12 
GPa. (b) Calculated lower bounds of correlation lengths for 
the superlattice. Separate points at a given P represent scans 
through different superlattice orders. All measurement tem- 
peratures were between 3.5 and 5.8K. 



are the atomic form factors at the superlattice and (3, 
1, 1) positions, respectively. Given the domain degener- 
acy in each superlattice order and uncertainty over the 
direction of S, we evaluate |q-£/2| 2 by averaging the di- 
rection of superlattice displacement 8 uniformly over the 
full 47T solid angle. Further averaging was carried out 
over the four (1/2, 1/2, l/2)-type superlattice domains. 
In this way we estimate S — 2.3 ± 0.7 10~ 3 A, which is 
a factor of 4 x 1CP 4 smaller than the lattice constant a. 
Using a value for the lattice force constant appropriate 
to pyrite-structured 3c? transition metal dichalcogenides 
(~ 1 N/cm = 6.25 eV/ A 2 , Ref. [23]), we estimate that 
the elastic energy associated with the superlattice dis- 
tortion is approximately 0.017 meV per Ni atom. This 
is at least two orders of magnitude smaller than the in- 



sulating gap and the magnetic exchange coupling, and 
we therefore consider the superlattice to be a minimally 
intrusive representation of the underlying magnetic or- 
der. We also point out that the superlattice distortion 
cannot account for the insulating behavior. For a wide 
band model appropriate to MS2 (W ^> E g ) both the lat- 
tice and electronic energies scale as 5 2 , and we need only 
confirm that net quadratic coefficient is positive. Using 
Eq. (3.54) of Ref. [5] and inserting appropriate values we 
find that the quadratic coefficient lies in the range 2.4 - 
3.1 eV/A . Rather than coupling to the formation of an 
energy gap, the superlattice is likely driven by the varia- 
tion in exchange constant J with Ni ion displacement!^. 

Using the superlattice reflections as a measure of the 
M2 order, we are able to track the magnetism and 
the crystal lattice through the pressure-driven insulator- 
metal transition. The only high-pressure magnetic scat- 
tering study published to date by Panissod et al. (Ref. 
[IS]) was limited to P < 2.9 GPa at 4.2 K with no dis- 
appearance of antiferromagnetism observed. This study 
identified the insulator-metal phase boundary with a sus- 
pected lattice discontinuity at 1.3 ± 0.4 GPa, leading to 
the conclusion that the magnetism is continuous across 
the insulator- metal transition 16 . Here we observe the 
(1/2, 1/2, l/2)-type superlattice distortion at every pres- 
sure from to 2.84 GPa (Fig. 2a) at base tempera- 
ture. Above 2.84 GPa the superlattice vanish es in all 
samples (Fig. 2c-d). Published account s ' 15 * 17 ! and our 
own transport measurements place the critical pressure 
for the insulator-metal phase boundary in the range 2.2 
- 3.1 GPa. This range brackets the upper limit of the ob- 
served superlattice diffraction. It is therefore natural to 
conclude that the disappearance of the low temperature 
M2 magnetic state, the structural phase transition (see 
below), and the insulator-metal transition all coincide at 
P ~ 2.9 GPa. 

The relationship between the correlated insulator and 
M2 antiferromagnetism is further revealed by considering 
the phase boundary T/V2 as a function of lattice constant 
for both NIS2 under pressure and NiS2-a;Se x in the P-x-T 
phase diagram (Fig. lc). Antiferromagnetic coupling of 
correlated electrons results from superexchange through 
the S ligand fields^. This coupling grows stronger as the 
lattice constant is reduced^, whether by chemical substi- 
tution or by applied pressure, continuing smoothly and 
continuously across the substitution-pressure interface. 
For the NiS2- 2: Se a ;(P, x) system, T/v 2 increases as the lat- 
tice constant shrinks, consistent with the superexchange 
interaction. The M2 magnetism and the associated su- 
perlattice distortion should be considered as byproducts 
of electron correlation and are not by themselves respon- 
sible for driving the insulating state. Important evidence 
for this also comes from electrical resistivity data which 
show that the Arrhenius activation energy is unchanged 
on cooling through T/V2 (Fig. lb). The insulating energy 
gap is thus well established before the formation of the 
M2 phase, reflecting the fact that the energy scales of 
the charge and magnetic interactions are well separated. 
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The Hubbard U is comparable to the Ni 3d e g bandwidth 
W = 2.1 eV, the magnetic exchange coupling is compa- 
rable to the transition temperature kgTpf 2 = 2.6 meV, 
and the supcrlatticc distortion energy scale is an even 
smaller 0.02 meV. 

Careful study of the lattice structure reveals additional 
information about the nature of the transition. High- 
resolution 9-26 longitudinal scans at all orders of the 
superlattice reflection are consistently close to resolution- 
limited (Fig. 2a), giving M2 correlation lengths greater 
than 1000 A(Fig. 2b). This is consistent for all pressures 
within the phase boundary. For the fee lattice peaks, 
longitudinal scans (Fig. 3) reveal a more complicated 
picture. All measured lineshapes are resolution-limited 
for < P < 0.85 GPa. However, beginning at 0.85 
GPa we observe multiple splitting of diffraction peaks in 
almost every sample at every pressure. This indicates 
the emergence of structural domains of reduced symme- 
try. Importantly, the superlattice reflections remain res- 
olution limited before disappearing entirely above 2.84 
GPa (Fig. 2a). The contrast between the multiply-split 
fee Bragg peaks and the sharp superlattice reflections is 
proof of phase coexistence between 0.85 GPa and 2.84 
GPa. We note that the onset of the high-pressure struc- 
tural phase that we observe at 0.85 GPa may explain 
the phase boundary at 1.3 ±0.4 GPa claimed in previous 
neutron scattering worM^. 



IV. LATTICE STRUCTURE UNDER 
PRESSURE 

The crystal symmetry of the high-pressure metallic 
phase can be determined cither with single crystal re- 
finement of a single-domain specimen or with powder re- 
finement of a polycrystalline sample with full knowledge 
of symmetry split peaks. We identified a single-domain 
sample in the high pressure phase at 2.9 GPa. Six diffrac- 
tion orders were measured and a least-squares refinement 
reveals an almost monoclinic structure with lattice pa- 
rameters (to a 95% confidence level) a — 5.5852(22) A, 
b = c = 5.6021(6) A, a = 89.984(8)°, /3 = 89.930(18)°, 
and 7 = 89.967(13)°. The values of (a-c)/a and (90°-/3) 
measured at P = 2.9 GPa are consistent with measure- 
ments at 3.86 GPa (Fig. 3), where the four-fold splitting 
of the (2, 2, 0) peak only can be explained by a symme- 
try of monoclinic or lower. Constraining the symmetry 
to monoclinic, we obtain (to a 95% confidence level) a 
= 5.5748(9) A, b = c = 5.5853(6) A, /? = 89.949(1)° 
at 3.86 GPa. The four-fold splitting of (2, 1, 1) at 0.85 
GPa (Fig. 3) is consistent with this picture, as it cannot 
be explained by a single phase of symmetry higher than 
monoclinic. 

Notably, symmetry reduction on passing into the 
metallic phase is the opposite of what is observed in many 
other transition metal oxides including the prototypical 
Mott-Hubbard system V2O3, which is a rhombohedral 
metal and a monoclinic insulator. For a non-interacting 
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FIG. 3: (color online), (a) Longitudinal 6 — 26 scans of the (2, 
1, 1) lattice reflection at 0.85±0.10 GPa, the lowest P at which 
phase separation was observed, (b) Rocking curve recorded 
at q = 2.721A for the same (2, 1, 1) reflection. The colored 
bars indicate the three different central 6 values used in the 
scans in panel a. (c) Longitudinal scan of the (2, 2, 0) lattice 
reflection at 3.86 GPa where only the high P phase is present. 
The four-fold splitting of (2, 2, 0) indicates a lattice symmetry 
of monoclinic or lower in the metal, (d) Lattice constants vs. 
pressure for pure NiS2 at T — 5 K. Both the low P and high P 
phases were observed in the phase separated region, but with 
insufficient detail to determine the lattice parameters a, b, c of 
the high P phase. The lattice constants ao of the low P phase 
were determined from the d-spacings of both the superlattice 
and cubic fee reflections. 



bandstructure, a red uction of lattice symmetry favors in- 
sulating behavioJ^El. The observation here for N1S2 that 
the symmetry is reduced in the metal therefore empha- 
sizes the role played by electron correlations in the insu- 
lator. 



V. CONCLUSION 

Our results address longstanding debates over the role 
of magne tism and crystal structure at the insulator- metal 
transitio n 11 * 2 * 4 * 5 * 24 -* while at the same time raising ques- 
tions about quantum phase transitions in the presence of 
strong electron correlations. The broad regime of phase 
coexistence that we observe while tuning U/W adds to 
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a growing list of correlated electron systems that exhibit 
phase coexistence around a first order quantum phase 
transition^. We have established that the magnetic, su- 
perexchange interaction cannot account for electron lo- 
calization in the insulator and that strong electron cor- 
relations drive the insulator-metal transition even in the 
presence of a structural distortion. Magnetotransport 
measurements in the compressed metal are required to 
probe the evolution and gapping of the Fermi surface, as 
well as the role of quantum fluctuations, as the transition 
is approached from above. 
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